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Abstract 

Wormlike micelle solutions are submitted to small-amplitude oscillatory shear 
superimposed to steady shear in the shear banding regime. By imposing a shear 
oscillation, the interface between high- and low-shear regions oscillates in time. A 
two-fluid semi-phenomenological model is proposed for superposition rheology in 
the shear banding regime, which allows us to extract a characteristic velocity for 
the interface dynamics from experiments involving only a standard rheometer. 
Estimates of the stress diffusion coefficient T> can also be inferred from such 
superposition experiments. The validity of our model is confirmed by directly 
recording the interface displacement using ultrasonic velocimetry. 
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Introduction 



During the last decade, wormlike micelle solutions have become a model 
system to study the so-called "shear banding" phenomenon. Depending 
on the concentration, most of these surfactant systems constituted of long, 
cylindrical, semi-flexible aggregates undergo a shear-induced transition from 
a state of entangled, weakly oriented micelles to a state of highly aligned 
micelles above some critical shear rate 77. Such a transition is strongly 
shear-thinning since the viscosity of the aligned state can be orders of mag- 
nitude smaller than the zero-shear viscosity of the system. Under simple 
shear and above 7/, the system spatially separates into coexisting bands 
of high and low viscosities corresponding respectively to the entangled and 
aligned states. As the shear rate is increased above 7/, the shear-induced 
structure progressively expands in the sample along the velocity gradient di- 
rection until the system is fully aligned at some shear rate 7Ar (in this work 
I and N respectively stand for isotropic and nematic in reference to the 
isotropic-to-nematic transition although the precise structure of the shear 
bands is still unclear). The rheological signature of shear banding is the 
existence of a horizontal plateau at a constant shear stress a = cxc in the 
shear stress vs shear rate constitutive curve (7(7), which extends from 7/ to 
7Ar. The present paper is restricted to the shear banding scenario described 
above and referred to as "gradient banding" in the literature. Another sit- 
uation known as "vorticity banding" may also occur in wormlike micelles, 
where the system separates into bands bearing different stresses stacked 
along the vorticity direction, corresponding to a vertical portion in the flow 
curve, i.e. to a shear-thickening transition. A recent revi ew of the sp ecific 



rheological properties of wormlike micelles is available in [Berret 2005 



The first experimental evidence for a stre ss plateau in nonlinear rhe 
ological measurements was provided by [Rehage and Hoffmann 1991 



on the CPCl-NaSal system. Further research effort established 
the general ity of this peculiar featu re o n other wormlike micelle 
systems (Berret and co-workers 19941 . iBerret and co-workers 1997 



ISoltero and co-workers 19991 . Theoretically, shear banding was first 
interpreted in th e framework of noneq uili brium phase transitions in 
liquid crystals (Gates and Milner 19891 . lOlmsted and Goldbart 1990l . 
lOlmsted and Goldbart 1992| . Specific features of the wormlike mi- 



celles such a s polymer- like behaviour and reversible breakage were then 
included by [Spenley and co-workers 1993| i n connection with the nori 
linear rheology of conventional polymers [Gates and co-workers 1993 



These two different approaches led to a theoretical debate about 
non-monotonic constitutive equations and shear banding seen ei- 
ther as a mechanical instability or as a nonequilibrium phase tran- 



sition [Schmitt and co-workers 1995 . Spenlev and co-workers 19961 . 

Olmsted and Lu 1997 . Porte and co-workers 1997i r Theoretical and 



numerical works later focused on including stress diffusion to ac- 
count for a unique s t ress selection and for the band dynamics 
Dhont I999I . Yuan I999I . lOlmsted and co-workers 2000| and on study ing 



the effects of flow-concentration coupling Fielding and Olmsted 2003a| or 
the possible instabilities inherent to the models Fielding 2005 . 

From the experimental point of view, phase separation under shear was 
ascertained for the first time by flow birefringence which showed the co- 



2 



existence of bands of weakly oriented and highly anisotropic material in 
sheared CTAB solutions close to an isot r opic-to-nematic equilibrium tran - 



sition [Cappelaere and co-workers 19951 . iMakhloufi and co-workers 1995 



Early nuclear magnetic resonance measurements confirmed the existence 
of inhomogeneous flows and the presence of differently sheared regions 



characterized by different order parameters Mair and Callaghan 19961 . 



Britton and Callaghan 19971 . [M air and Callaghan 1997l | but it is not un- 
til recently that the simple shear banding scenario described above re- 
ceived full experimental validation from light scattering and particle track- 
ing velocimetry in the CPCl-NaSal system Salmon and co-workers 20031 . 



Mendez- Sanchez and co-workers 20031 . iHu and Lips 2005| . In particular the 



so-called "lever rule" which, in strong analogy with first-order equilibrium 
phase transitions, gives the proportion a of the aligned state as a function 
of the shear rate 7 along the stress plateau: 



7 



«)7/ + Ol^N, 



(1) 



appears as a rather robust f e ature provided that steady s tate is reached 



Salmon and co-workers 20031 . iLerouge and co-workers 2004 



Thus, 
debate, 

established [Lopez- Gonzalez and co-workers 2006 



although the exact nature of shear bands is still under 
the coexistence of differently sheared bands is now well 

Most of latest 



work on shear bandi ng has concentrated on the l ocal flow dynamic s 



during transients 
on velocity and birefringence 
t he steady scena rio described 



Lerouge and co-workers 2004 iHu and Lips 2005 



fluctuations 
above 



Becu and co-workers 2004. 



and departures frorn 



Holmes and co-workers 2003 



Lopez-Gonzalez and co-workers 2004 



Lee and coworkers 2005. Yesilata and coworkers 2006 , on interface stabil 



ity [Lerouge and co-workers 2006 



dynamics 



and on m odelling such spatio-temporal 
Radulescu and co-workers 2003 . Fielding and Olmsted 2003bl . 



Fielding and Olmsted 2004 iFielding and Olmsted 2006 . 

In this pap er we propos e to use the parallel superposition technique 
introduced by [Booij 1966al | to investigate shear banding in wormlike mi- 
celles and more precisely to access the dynamics of the interface be- 
tween shear bands. In our opinion, the interest of superposition rheol- 
ogy has been overlooked in the literature. In particular, only a very lim- 
ited number of papers are devoted to superposition measurements in com- 
plex fluid s that show strong fiow-microstruct ure coupling, e.g. associative 

{ )olymers Tirtaatmadia and co-workers 199"7j | or liquid crystalline polymers 
Grizzuti and Maffettone 2003| . Our aim is to show how this technique, 
which is available on most rheometers, can be used to access the dynamical 
behaviour of shear bands, without having to rely on involved techniques 
as described above. We first recall the principle of superposition rheology 
and illustrate it in the case of wormlike micelles sheared below 7/, i.e. in 
the homogeneous, entangled state. Then a two-fluid semi-phenomenological 
model is described for superposition rheology in the shear banding regime in 
the simple case of infinite parallel plates. This model is extended to account 
for experimental geometries, namely cone-and-plate, Couette, and Mooney- 
Couette geometries. The corresponding calculations are gathered in the 
appendix. Finally our model is probed experimentally on the well-studied 
wormlike micellar system CPCl-NaSal through superposition rheology and 
compared to direct measurements of the interface dynamics using ultrasonic 
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velocimetry in Couette geometry. The results are further discussed and in- 
terpreted in terms of the stress diffusion coefficient a key parameter in 
recent theoretical approaches of shear banding. 



ONE-FLUID SUPERPOSITION RHEOLOGY 



Superposition rheology as first introduced by jBooij 1966a| is the addi 



tion of a small-amplitude oscillatory shear to a main steady shear. The 
oscillatory shear can be either parallel or perpendicular to the steady shear. 
Superposition allows one to probe the dynamical response of a shear-driven 
system and to generalize the notions of viscoelastic moduli to far-from- 
equilibrium conditions through a perturbation analysis. The properties of 
the "superposition moduli" and the relations hips between "parallel moduli" 
and "orthogonal moduli" were discussed by [Vermant and co-workers 19981 



and [Phont and Wagner 200 1| and applied to polymer solutions and col 



loidal suspensions, respectively. Here we focus on parallel superposition 
which is now available as an option on most recent commercial rheometers. 



A. Notations for one-fluid superposition rheology 

Let us first introduce the various notations for superposition rheology. 
In the following we shall use complex notations and assume that the shear 
rate reads 

7 = 7i+72e'"*- (2) 

Both 7i and 72 are taken to be real and positive. If 72 corresponds to a 
perturbation to the steady shear in the linear regime, the shear stress can 
be written as 

a = (Ti + ^26^"* , (3) 

where ai is real and (T2 is the complex amplitude of the oscillatory part of the 
shear stress. The issue of specifying which variable is controlled and which is 
measured will be addressed below and further discussed in sect . IIIBI From 
eqs. (E]) and ([3]), two apparent viscosities are defined as 

V = 4 

71 

"11 = ^- (5) 

Since the oscillatory part of the shear is only a linear perturbation of the 
steady component, rj depends only on 71 and reduces to the standard shear 
viscosity found when 72 = 0. On the other hand the complex viscosity 
?7il depends on both 71 and u, and allows one to explore the dynamical 
behaviour of the shear-driven system. 

When both 71 and 72 tend to zero, one should recover the usual complex 
viscosity ri*{Lj) so that 

lim r7j|(^,7i) = r/*(^) . (6) 

7i->0 " 
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FIG. 1: Linear rheology of an 8% wt. CPCl-NaSal solution: storage modulus G' 
(•) and loss modulus G" (o) versus frequency lo. The solid lines correspond to 
an Oldroyd-B fluid (eq. ^) with 770 = 122 Pa.s, n = 0.87 s, and T2 = 0.60 ms. 

Another useful limit is found by considering vanishing frequencies for a 
finite 71. In that case, 72 and (J2 become steady perturbations so that 
eq. (IS]) reduces to r/jj = dcri/d7i, which leads to 

d?7 

lina r/|[(w,7i) = r/(7i) +71 — (71) . (7) 
w^o " d7i 

B. Conventional rheology of wormlike micelles in the low-shear 
regime 

As already reported many times in the literature, semi-dilute so- 
lutions of wormlike micelles present an almost perfect Maxwellian be- 
haviour in the linear regime [Rehage and Hoffmann 1988 |. Such a strik- 



ing feature was predic ted and explained in terms of a reaction-diffusion 
model by [Gates 19"87 |. However, at high frequencies, significant devia- 



tions from the Maxwell rnodel may occur due to fas t relaxation modes 
Fischer and Rehage 199?! . lYesilata and coworkers 2006| . Thus a more thor- 



ough description of the low-she ar rheology of polymer-like rn icelles is pro- 
vided by the Oldroyd-B model Oldroyd 19531 . lOldroyd 1955l | whose linear 



complex viscosity 77* (tu) and nonlinear shear viscosity 77(7) read 

V ^ = ^0 T— , (8) 

1 + lUTi 
1 + (S27)' 

^(7) = ^0 , (9) 

where ri, T2, Si, and S2 are characteristic times. 

In the following, we focus on a wormlike micelle solution made of 
cetylpyridinium chloride (CPCl, from Aldrich) and sodium salicylate 
(NaSal, from Acros Organics) dissolved in brine (0.5 M NaCl) with 




FIG. 2: Nonlinear rheology of an 8% wt. CPCl-NaSal solution: shear stress 
a versus shear rate 7 (a) in linear scales and (b) in logarithmic scales (right). 
The solid line corresponds to an Oldroyd-B fluid (eq. ([9])) with 770 = 122 Pa.s, 
si = 0.59 s, and S2 = 0.13 s. The dashed line is the best fit of the high-shear 
branch by a Bingham fluid a = o'B + 'ilB'y with aB = 91.7 Pa and i]b = 1.13 Pa s. 
The shear banding regime extends from 7/ ~ 2.2 it 0.2 s~^ to 7Ar — 7.4 it 0.4 s~^. 

a fixed concentration ratio [NaSal]/[CPCl]=0.5 and a total surfactant 
concentration of 8% wt. (unless stated differently) as described by 
(Rehage and Hoffmann 19881 . iBerret and co-workers 1997 ] . The working 



temperature is T = 21°C. Figure [T] shows the linear viscoelastic moduli of 
our micellar solution measured in the Mooney-Couette geometry described 
below (see sect. IIII A|) with a standard stress-controlled rheometer (ARIOOO, 



TA Instruments). All the experiments in the present work were performed 
under controlled shear stress. Both G' and G" are very well described by 
the Oldroyd-B model (eq. ([8]) with G' + iG" = iurj*) which captures the 
departure of G" from the scaling at high frequencies. 

The constitutive curve a vs 7 of the same micellar solution is shown in 
fig. [2j As expected the fluid is weakly shear-thinning below 77 ^ 2.2 s~^. 
Above 7/ very strong shear-thinning is observed and the stress saturates 
at a plateau value — 100 Pa. This corresponds to the shear banding 
transition. The solid line in fig. [2] shows that the nonlinear rheological 
behaviour of our fluid in the low-shear regime is rather well captured by the 
Oldroyd-B model (eq. Q). 



C. Superposition rheology of wormlike micelles in the low-shear 
regime 

The super position rheology of an Oldroyd-B fluid was computed by 
|Booij 1966bl | and leads to 

^11 (^^71) _ 1 - nT2Uj\l + sj'jl) + {3sl ~sl + sjsl'jD'yf + iuj [n + + [nsl + rgg^) 7?] 
Vo (1 + s?7f) [(1 + iunf + sj^f] 

(10) 
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FIG. 3: Superposition rheology of an 8% wt. CPCl-NaSal solution in the low- 
shear regime: | 7yjj(a;,7i) | versus uj for (a) 71 = 0, (b) 0.025, (c) 1.09, and 
(d) 1.49 s~^. The solid lines correspond to an Oldroyd-B fluid (eq. (jlOp ) with 
r/o = 122 Pa.s, si = 0.59 s, S2 = 0.13 s, ti = 0.87 s, and T2 = 0.60 ms. 

It is easily checked that eqs. ([6]) and ([7]) are recovered from eq. (ITOl) when 
the limits 71 — and u; — > are considered. 

Figure [3] presents superposition measurements in the low-shear regime. 
The steady-state shear rates 71 indicated in the caption of fig. [3] (and later 
figs, m El and [T3l) are the values measured by the rheometer. In all our 
experiments, the amplitude of the oscillatory part of the shear stress is 
fixed to (J2 = 0.5 Pa, except for sect. IIIIBI where (T2 = 1 Pa. Figure [3] 
clearly shows that the four parameters inferred from figs. [U and [2] yield a 
good description of ?7*(u;,7i) for all 71 < 1.5 s"'^ when used in eq. ( ITOi) . 

II. TWO-FLUID SUPERPOSITION RHEOLOGY: THEORETI- 
CAL PREDICTIONS 



The above results obtained in the low-shear regime prompt 
us to use the superposition technique in the shear banding 
regime. Indeed the CPCl-NaSal system is known to separate into 
weakly and highly sheare d bands as described in t he introduction 
iBerret and co-workers 1997. Porte and co-workers 1997 . Previous 



work has shown that the shear banding phenomenon is rather simple 
in this particular system: the proportion of shear-induced struc- 
ture is given by the lever rule ([I]) and no wall slip is detected 
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FIG. 4: Superposition rheology of an 8% wt. CPCl-NaSal solution in the shear 
banding regime: | ?/||(u;,7i) | versus w for (a) 71 = 2.5 and (b) 7.14 s~^. The solid 
lines correspond to an Oldroyd-B fluid (eq. (jlO|) ) with r^o = 122 Pa, s\ = 0.59 s, 
S2 = 0.13 s, Ti = 0.87 s, and T2 = 0.60 ms. 



Salmon and co-workers 20031 . iHu and Lips 20051 



8% wt. CPCl-NaSal solution, 
7 > 7/ ~ 2.2 s~^ and a = cr^ 

shear stress is in good agreement w ith the prediction ar 

Go = Vo/'^i the plateau modulus [Spenley and co-workers 1993 



In the case of our 
the shear banding transition occurs for 
~ 100 Pa. The value for the critical 

0.67Go, with 
Note 



however that the stress plateau is not perfectly fiat at in fig- El This is 
most probably due to the curvature of the Mooney-Couette geometry which 
induces a significant slope of the consti tuti ve curve in the shear bandin g 
regime [Radulescu and Olmsted 2000l . ISalmon and co-workers 2d03 |. 
Such a s lope may also arise from fiow-concentration cou 
pling 



(Schmitt and cq-workers 1995 



Olmsted and Lu 1997 



Fielding and Olmsted 2003a| . However in the absence of clear exper- 
imental evidence for such a mechanism in the literature on the system 
under study, we shall not refer to concentration coupling effects thereafter. 
In any case, the slope in the fiow curve makes it hard to distinguish between 
the "stress plateau" and the homogeneous high-shear regime where the 
system is fully aligned (see also appendix IA2[) . 

Figure |4] shows that the Oldroyd-B model used in the low-shear regime 
completely fails in describing the complex viscosity 77^(0;, 71) when 71 > 7/. 
This is a strong indication that the system enters the shear banding regime 
and that the model for superposition needs to be modified. In the follow- 
ing, we discuss a two-fiuid model for superposition rheology in the presence 
of shear banding. This simple model is presented for various geometries, 
from the most simple geometry (infinite parallel plates) to the more com- 
plicated one actually used in our experiments (Mooney-Couette geometry). 
The detailed calculations for experimental geometries are presented in the 
appendix. 
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A. Infinite parallel plates 



Let us first consider the case of two unbounded parallel plates in transla- 
tion separated by a gap e. We assume that the fluid separates into bands of 
"isotropic" (I) and "nematic" (N) material. In a superposition experiment 
in the shear banding regime, the steady component of the shear stress is 
fixed to (Ti = Gc and the steady component of the shear rate in the isotropic 
(resp. nematic) material is simply 7/1 = 77 (resp. 7Ari = 7Ar), where 77- 
and 7Ar are the limits of the stress plateau. This leads us to generalize the 
notations introduced in sect. II Al to the two-fluid case: 







(Tc + (726 , 


(11) 


at) 




7l + 726 , 


(12) 






ll + 7/26 , 


(13) 


7iv(t) 




In + 7iV2e*'^* , 


(14) 


Vi = — 

7/ 


and 


f^2 

^11/ = — , 
" 7/2 


(15) 


Vn = — 

In 


and 


0-2 


(16) 



Since the steady shear in each phase is fixed to 7/ and 7Ar respectively, r]i 
and ?77v are two constants and ri*^j and r^jj^ depend only on uj. In other 
words vi and vn are the apparent viscosities of the isotropic and nematic 
materials under a steady shear stress ai = Oc-, while r(^j and rf^^ correspond 
to the dynamical behaviours of the two phases for a\ = (Tc- 

In the absence of wall slip, the lever rule ([1]) is a mere conse- 
quence of the continuity of t he velocity at the interfaces between bands 
Salmon and co-workers 2d03l | and the proportion of shear-induced struc- 



ture a{t) obeys 

a{t) = ai + ^26^'^* , (17) 
7(t) = (1 - «(t))7/(t) + a(t)77v(t) , (18) 
7i = (1 -01)77 + ai77v. (19) 

Note that the fact that the "instantaneous" lever rule (fT8l) applies at all 
times actually results from a steady-state approximation of the Navier- 
Stokes equations, i.e. from assuming that pdtv -C dzCr. This assumption 
will be checked below a posteriori. Moreover eq. (fT9ll leads to the lever rule 
for the apparent viscosity 

1 = ^ = 1^ + ^. (20) 
?7 Vi Vn 

Using the above notations and restricting the analysis to linear response, 
it is easily shown that 

(T2 cr2 / 1 1 \ 

72 = (1 - ai) — + ai—— + a20-c • (21) 

V\\i V\\n \Vn VI J 

In order to get an expression for r/jj = (72/72 in the shear banding regime, 
we need to link a2 and (T2. We chose to use the reaction-diffusion model 
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proposed by [Radulescu and co-workers 1999l | that assumes the existence 



of a single band and shows that the interface between the isotropic and 
nematic regions moves at a velocity c that only depends on the difference 
(T(rc) — CTc, where Vc is the position of the band, and vanishes for cr{rc) = (Tc- 
More precisely, if one assumes the shear-induced structure to be located 
from r = to r = rc{t) = a{t)e, with r being the coordinate across the gap, 
the model predicts 



da a(r) — Or dc 



dt e da 



a{r) -a^CQ _ a2e"^^ Cq ^^2) 



(Tc e (Tc 



where we have introduced the characteristic velocity cq defined by co/o"c = 
dc/do"!^^. Equation (fTTIl then leads to 

0"2 Co 

^2 = — -. — • (23) 
(Tc luje 

Experimentally the amplitude of the oscillations of the interface position 
rc{t) should be accessible through the time-resolved velocimetry techniques 
mentioned in the introduction provided that the spatial resolution is fine 
enough. An estimate of Cq from direct measurements of a2 using ultrasonic 
velocimetry will be presented in sect. IIII B[ Finally inserting eq. (!23l) into 
eq. (12T!) and using the definition ([5]) yields 

]_ ^ 1 - ai _^ _^ 7jv - ii Co 
Vii Viii ViiN iuje' 

Equation (1241) shows that the complex viscosity in parallel superposition 
involves two terms: 

- = ^ + ^, (25) 
1 7^-7/ CO ^26) 



The first term 77^(71, tu) corresponds to the "steady" lever rule fl20l) applied 
to the complex viscosities r/yj- and r/jj^ and depends on both 71 (through ai) 

and ui (through 77^^- and ?7||^). In principle rj'^^jiuj) and 77^^(0;) are accessible 
through superposition measurements in the homogeneous states at 71 = 7/ 
and 7i = 7Ar respectively (or at least by extrapolation of ?7ij(7i,ci;) when 

7i — > 77 and 71 — > 7]^), so that ?7i(ciJ,7i) is known once ai is known via 
eq. (USD. 

The second term ri£,^{uj) accounts for the dynamics of the interface be- 
tween the two bands and does not depend on 71. Note however that this 
"dynamical" term depends on the geometry since e shows up in eq. fl26l) . 
Since good approximations of 7/, 77V, and (Tc are given by nonlinear rhe- 
ological measurements, the only unknown in eq. fl26|) is Cq. We conclude 
that superposition rheology in the shear banding regime should provide an 
experimental means of probing the dynamics of the interface between shear 
bands through the measurement of Cq. In practice the various parameters 
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involved in eq. (12^ are not that easy to extract from independent measure- 
ments. As already pointed out the limits of the stress plateau are not always 
clear (see fig. [2]). But the largest difficulty probably lies in getting a good 
approximation for the dynamical behaviour of the shear-induced structure 
?7j|jy(ci;) from superposition measurements at 71 > 'Jn- Indeed the high-shear 
branch of the flow curve is sometimes impossible to access due to flow in- 
stabilities t hat tend to expel the sample from the measuri ng tool at high 
shear rates Berret and co-workers 1997I . IHu and Lips 2005 |. Still one may 



argue that 77*^(0;) could also be inferred from superposition experiments in 
the shear banding regime by looking at the dependence of ?7*(ci;,7i) on ai 

in eq. fl2^ . We shall further discuss this point below in sect. IIII CI 

A simple way to overcome the difficulty raised by ?7i|^ is to focus on the 

limit tti 0, i.e. just at the onset of shear banding. In this limit eq. ( l24l) 
becomes 

liml = ^ + ^:^^, (27) 

where Cq is the only unknown parameter since r/y^- is known from measure- 
ments in the low-shear regime. Experimentally, ai is varied for a given 
frequency ut. The value at the origin of the linear regression of l/?7* vs ai 
is then eq. ( 1271) . from which l/rjo^oi^) determined. Finally a linear fit of 
rjDooi'^) '^^ ^ yields cq. This fitting procedure will be tested in sect. IIII AI 

To conclude this discussion of eq. (l2ll) . let us check the validity of the 
"instantaneous" lever rule ( ITSl) . In the case of oscillating velocity and stress 
fields, neglecting the time derivative in the Navier-Stokes equation is equiv- 
alent to setting puJV2 -C o"2/e where V2 — 726. In the low-frequency limit, 
eq. IHM yields 77* = (T2/72 — iujead cq{^n — ii) so that the steady-state 
approximation holds if 

pcoe < 1 . (28) 

With the typical values p = 10^ kg.m"'^, cq < 1 mm.s~^ (as will be checked 
experimentally below), e = 1 mm, 7Ar — 7/ = 10 s~^, and cXc = 100 Pa, the 
left-hand side of eq. fl28|) is about 10~^. Hence the approximation holds at 
least in the low- frequency limit of eq. fl2^ . More generally, the steady-state 
approximation reads pue"^ <^| ?7j| |. It can be checked from figs. [3] and Hit hat, 

for the highest frequencies achieved in our experiments {u ~ 100 rad.s"^), 
one always keeps | ?7|| |^ 1 Pa.s, so that, with pue'^ < 0.1, the approximation 
remains valid. 



B. Experimental geometries 

Standard experiments use cone-and-plate or Couette geometries (or their 
combination known as the Mooney- Couette geometry). The changes that 
the use of such geometries induces in eqs. f l2^ - fl26l) are described in details 
in the appendix. It is shown in sect. I A II and I A 21 that the expressions found 
for rj* in both the cone-and-plate and the Couette geometries can be written 
in forms similar to eq. (12^ . 
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In particular, in both cone-and-plate and infinite parallel plates, it is 
seen from eq. (123|) that 02 diverges at low frequencies for fixed (T2, i.e. 
under controlled stress. Such a behaviour is a direct consequence of the 
highly nonlinear fluid response under controlled stress in fiat geometries, 
where jumps between the two shear branches of the flow curve are expected. 
Therefore, to ensure that the experiments are conducted inside the stress 
plateau for all frequencies, superposition rheology in the cone-and-plate 
geometry requires to work under controlled shear rate, so that 0:2 (and thus 
72 through eqs. fll7p - ffT^ ) always remains a linear perturbation of the steady 
shear. Since a controUed-stress rheometer is used in the present work and 
since ultrasonic velocimetry is not available in the cone-and-plate geometry, 
we shall rather focus on the Couette geometry where the divergence of 02 
does not occur. 

More precisely, sect. lA 21 shows that in the "small-gap approximation," 
i.e. when the gap e is small enough compared to the radius Rq of the 
inner cylinder, the case of a Couette geometry reduces exactly to the case 
of infinite parallel plates provided that riD^ is replaced by 



1 7Af - 7/ Co 



(29) 



This corresponds to the zero-order version of eq. (1A14I) . i.e. it assumes that 
both the term of order e/ Rq in eq. ( 1A16I) and the first-order corrective term 
1/rjQ given by eq. ( 1A17P can be neglected. As seen on eq. ( lAllI) the curvature 
of the Couette geometry prevents 0:2 from diverging at low frequencies so 
that superposition measurements can be performed under controlled stress. 

Finally sect. IA3I shows that the case of a Mooney- Couette geometry of 
height h can be handled by considering the proportions = (1 + -Ro/2/i)~^ 
and = 1 — eco of the surface respectively covered by the Couette (co) and 
by the cone-and-plate (cp) geometries relative to the total surface. In par- 
ticular, ?7* is given by the following average of the corresponding viscosities 

^\\co and r/*^p : 

V\\= ^coV\\co + ^cpV\\cp- (30) 

However to close the problem one has to specify the values of ai in the 
two parts of the geometry. As shown in the appendix, this leads to serious 
complications and the interpretation of superposition measurements in the 
Mooney-Couette geometry requires in principle the full knowledge of the 
dynamical behaviours r/*^ and r/jj^ of the high- and low-viscosity materials. 

To keep things analytically tractable and although this may be a crude 
approximation of the actual behaviour in the Mooney-Couette geometry, 
we shall assume that 

A = ^ + ^' (31) 

'^ll ^\\co ^\\cp 

which is consistent with eq. (!30l) only for e/i?o <^ 1- Using the effective cti 
found in sect . lA 31 and given by 



Ro ( hi-ecp(Jc A 71-7/ .ooN 
tti = Ceo — \\ 1 +ecp^ (32) 



e \ V ^co<7c J In -li 
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FIG. 5: (a) Real and (b) imaginary parts of l/?7j| (71 , w) versus ai deduced from 
eq. (|A2ip (o), from eq. ()A22p (v), and from eq. ()A23P (♦). The solid lines are the 
best linear fits of the ♦ data while the dotted lines show the linear fits obtained 
using the o and V data. The frequency is (J = 0.26 rad.s"^. The fluid under 
study is an 8% wt. CPCl-NaSal solution. 

together with eq. flA14p at zero order in e/i?o for r/y^^ and eq. flA6p with 
e = i?o tan /3 for 77*^^, one finds 

Equation fl33|) is exactly eq. fl2ll) up to a corrective frequency-dependent 
term on cq that accounts for the Mooney-Couette geometry. Thus under 
the above assumptions we may still use the data analysis procedure de- 
scribed above in sect. Ill AI (see eq. fl27|) ) on experimental data recorded in 
the Mooney-Couette geometry. 

III. TWO-FLUID SUPERPOSITION RHEOLOGY: EXPERIMEN- 
TAL RESULTS 

A. Superposition experiments in the Mooney-Couette geometry 

Superposition experiments were performed in the shear banding regime 
on the previous 8% wt. CPCl-NaSal solution under controlled stress in 
a Mooney-Couette geometry with inner radius Rq = 24 mm, outer radius 
Ri = 25 mm, and height h = 30 mm. Using the notations defined above, 
this corresponds to — 0.7 so that we can not neglect the presence of the 
cone. The small-gap approximation holds since e/i?o — 0.04. A solvent trap 
is used to prevent evaporation and we checked that no significant change of 
the rheological properties of our micellar solution occurs over the ~ 8 hour 
maximal duration of our experiments. 
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FIG. 6: (a) | ri*^j{i^) \ measured closest to the onset of shear banding for 71 = 
1.7 s~^ < 7/ (o) and inferred from the fitting procedure based on eq. ([33]) (•). 
(b) I r]*^j^{uj) I measured in the high-shear state for 71 = 11.1 s~^ (o) and deduced 
from eq. ([33]) (•). The fluid under study is an 8% wt. CPCl-NaSal solution. 



Figure [5] shows the experimental l/?7y(ai) data obtained when varying 

the imposed steady shear stress ai {i.e. the average proportion ai of ori- 
ented phase) for a given frequency u. l/?7j|(ai) is inferred from the raw 
data ?7y(7i,Ci;) at fixed uj (see figs. HI and [6] for examples of such raw data). 
To test the robustness of the linear behaviour of l/ryjj vs ai expected from 
eq. f[33]) . the data were plotted against ai computed from eq. flA2ip alone 
{i.e. taking eco = 1 and neglecting the cone-and-plate part of the geometry, 
see o symbols), from eq. (IA22[) alone {i.e. taking ecp = 1 and neglecting the 
Couette part of the geometry, see V symbols), and from the full eq. flA23P 



with eco = 0.7 and = 0.3 (see ♦ symbols). The quality of the three 
linear fits are similar and the values of the slopes as well as the intercepts 
at ai = are all very close. We conclude that the linear behaviour pre- 
dicted by eq. f[33]) is indeed observed and that the way ai is computed is 
not critical. 

In order to use the extrapolation procedure proposed in sect. Ill AI for 
eq. ([27]) . ri*^j is taken to be the experimental value for the homogeneous 
fluid obtained closest to the onset of shear banding. The corresponding 
data are shown in fig. [6|(a) (o symbols, see also the discussion in sect. IIII C]) . 
We then calculate l/fj = l/?7j|(ai ^ 0) — foi' various frequencies uj 

ranging from 0.07 to 70 rad.s"^. The real and imaginary parts of fj are 
plotted as a function of uj in fig. [3 
If eq. ([33]) holds, one expects 



fl{uj) 



CTc 



luje 



In - li Co 
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LuRo 



+ 2e 



cp 



(34) 



If one further assumes that 2cq/ujRo <^ 1, then one should find a range of 
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FIG. 7: '^{fj) (•) and ^{fj) (o) versus uj in (a) logarithmic scales and (b) semi- 
logarithmic scales. The solid line is the best linear fit of 9(r/) by eq. (j35p with 



7/ = 2.2 s' 



7.4 s ^, cJc = 100 Pa, e = 1 mm, and cq = 0.1 mm.s ^. The 



fluid under study is an 8% wt. CPCl-NaSal solution. 



u for which Q{ri) ^ ^(v) 

$5(77) ^ 
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(35) 



cp 



Figure [Tl^a) shows that ^^(r]) ^ ^{v) for ~ 0.07-3 rad.s"^ in the ex- 
periment. The best hnear fit of '^{fj) vs u over this range of frequencies 
yields Co(l + ecp) = 0.13 ± 0.05 mm.s~^ so that Cq = 0.1 ± 0.04 mm.s~^. 
The large uncertainty 40%) on the determination of cq is mainly due 
to the uncertainty on 7/ and 77V and therefore on the calculation of ai. 
Since 2cq/ujRq ^ 0.003-0.1 for oj ^ 0.07-3 rad.s"^, the approximation lead- 
ing to eq. ( l35i) is justified a posteriori. These results were obtained with 
7/ = 2.2 s~^, 7Ar = 7.4 s~^, and = 100 Pa, which were estimated inde- 
pendently from nonlinear rheology as explained in sect. lA 21 However, at 
"high" frequencies (cu > 1 rad.s"'^), the terms induced by the curvature of 
the Mooney-Couette geometry are no longer negligible, so that first-order 
terms in e/i?o should be taken into account in eq. (1A14I) . This most prob- 
ably explains the observation of negative data for ^{fj) in fig. [7](b). 

In the last section of this paper we use ultrasonic velocimetry to directly 
access the dynamics of the interface during superposition experiments and 
check the validity of the above findings. These experiments were performed 
on a 6% wt. CPCl-NaSal solution (due to technical limitations involving the 
velocimetry setup and the 8% wt. sample). To allow for a direct comparison 
with velocimetry experiments, fig. [H] presents the analysis of superposition 
rheology measurements performed on the 6% wt. solution. The results are 
qualitatively the same as those for the 8% wt. sample shown in fig. [71 The 
estimate for cq in the 6% wt. sample is Cq = 0.31 ± 0.15 mm.s~^. 
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FIG. 8: Q{fj) (•) and ?R.{fj) (o) versus u> in (a) logarithmic scales and (b) semi- 
logarithmic scales. The solid line is the best linear fit of ^(r/) by eq. (|35p with 
7/ = 4.0 s~^, = 6.3 s~^, Uc = 68 Pa, e = 1 mm, and cq = 0.31 mm.s~^. The 
fluid under study is a 6% wt. CPCl-NaSal solution. 

B. Ultrasonic velocimetry during superposition experiments 



1. Velocity profile measurements 

Superposition experiments in the shear banding regime have shown the 
possibility of characterizing the dynamics of the interface between shear 
bands using only a standard rheometer. In this section, the above re- 
sults and model are confirmed using time-resolved local velocity measure- 
ments. To access the velocity field we used the ultras onic velocimetry 
technique described in [Manneville and co-workers 2004| . As shown by 



[Becu and co-workers 2004 1 this technique allows one to measure the ve- 
locity profile of shear-banding wormlike micelles in the gap of a Couette cell 
with a temporal resolution of about 1 s and a spatial resolution of about 
40 yum. 

Figure [9] shows a typical velocity profile v{r) measured in a 6% wt. CPCl- 
NaSal solutio n, where r is the distance from th e inner rotating cylinder. As 
explained in Manneville and co-workers 2004| . the fiuid was seeded with 



1% wt. hollow glass spheres (Sphericel, Potters Industries) of mean ra- 
dius 11.7 /im and density 1.1 in order to provide acoustical scattering. We 
checked that both linear and nonlinear rheological properties were not sig- 
nificantly affected by the addition of such acoustic contrast agents. Since 
the velocity profiles are recorded in the Couette part of the Mooney-Couette 
cell, we shall focus on the model developed in sect. IA2I Let us only recall 
here eq. ( lAllI) which gives the complex amplitude r2 = of the interface 
displacement: 

In the following, the steady shear stress is fixed to ai such that ai ^ 0.5. 
The steady-state velocity profile of fig. [9] clearly shows two linear parts that 
separate the gap into two shear bands of equal width where the apparent 
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FIG. 9: Velocity profile v{r) (o) recorded in a 6% wt. CPCl-NaSal solution 
at steady-state for cji = 70.5 Pa (which corresponds to ai ~ 0.5). The solid 
lines represent linear fits of the velocity profile in the high- and low-shear bands. 
Their intersection yields the position of the interface. The dotted line shows 
the velocity profile for a Newtonian fluid. 

viscosities differ by a factor of about two. Linear fits in the two shear bands 
yield the interface position ri = aic ~ 0.5 mm. 



2. Measurement of cq in a transient experiment 



In the framework of the model proposed by 
Radulescu and co-workers 1999j | the characteristic velocity cq can be 



deduced from transient velocity profile measurements. Indeed by suddenly 
decreasing the shear stress from a = cxi + (T2 to cr = cti at time t = and by 
measuring the evolution of the velocity profiles u(r, t) in time, we can easily 
track the interface position rc(t). Experimentally (J2 is fixed such that, 
by using eq. ( ]A9[) . r2(u; = = a2{uj = 0)e ^ 0.2 mm. In the small-gap 
approximation, the equation for the interface position reads: 

1 drc (T(rc) -ac 2 

— (n - r^it)) , (37) 



Co dt (Tc Ro 



which leads to 



r,(t) = ri + rs e'^'^"*/^" , (38) 

where r2(a; = 0) was simply noted r2. As seen in fig. [T0|, the position of the 
interface ^^(t) is well fitted by eq. fl38l) which yields Cq = 0.28±0.03 mm.s~^. 
Comparing with the results of the superposition measurements shown in 
fig. [HI one finds that both values are in quantitative agreement, which con- 
firms the relevance and the ability of superposition rheology to extract dy- 
namical information in the shear banding regime. Of course the uncertainty 
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FIG. 10: Position of the interface rc{t) versus time as the external shear stress is 
reduced from ai + cr2 = 71.5 Pa to cji = 70.5 Pa at t = 0. The solid line is the 
best fit by eq. (j38p with r2 = 0.17 mm and cq = 0.28 mm.s~^. The fluid under 
study is a 6% wt. CPCl-NaSal solution. 



on Co given by time-resolved velocimetry is much less than that of the su- 
perposition method (but at the cost of using a more involved technique and 
processing a large amount of ultrasonic data). 

The present analysis of the velocity measurements also neglects the 
first two stages of the band dynamics during the transient, namely low- 
shear band destabilization and interface reconstruction, as evidenced by 
|Radulescu and co-workers 20031 ]. These initial stages were shown to occur 
in typically 2 s which is of the order of the temporal resolution of our ve- 
locimetry experimen ts. Thus we only focus on the last dynamical step called 
"interface travel" in [Radulescu and co-workers 20031 ]. In particular, due to 
the existence of two early relaxation stages, one may argue that the initial 
position is ill-defined and that the interface position after reconstruction 
may significantly differ from ri +r2. This is the reason why r2 was actually 
left as a free parameter in eq. ( l38l) . 



3. Validation of the model for superposition experiments 

Now turning to the superposition experiment, we add an oscillatory shear 
stress of amplitude o"2 to a steady shear stress cxi, and follow the position of 
the interface in time for various frequencies. Figure fTTl presents the measure- 
ments of the local shear rates in the two bands and of the interface position 
rc{t) versus time for two different frequencies. As expected these various 
quantities oscillate in time and, in spite of some experimental scatter, fitting 
rc{t) by sine functions for various frequencies yields a good estimate for the 
amplitude r2{uj) = a2{uj)e. 

The dots (•) in Fig. [12] show the amplitude ] r2 ] of the interface os- 
cillations inferred from ultrasonic velocimetry for four different frequen- 
cies, while the solid line is calculated using eq. fl36l) with the value cq = 
0.28 mm.s~^ obtained from the transient experiment. The quantitative 
agreement between the experimental data and the calculated prediction con- 
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FIG. 11: (a) Local shear rates versus time in the nematic band (top) and isotropic 
band (bottom) along with the global shear rate recorded by the rheometer (mid- 
dle) for to = 0.05 rad.s"^ and (b) uj = 0.1 rad.s"^. Position of the interface rdt) 
versus time for (c) u; = 0.05 rad.s"^ and (d) to = 0.1 rad.s""*^. The solid lines are 
the best fits by sine functions. The fluid under study is a 6% wt. CPCl-NaSal 
solution submitted to stress oscillations of amplitude £72 = 1 Pa around the mean 
value o"! = 70.5 Pa. 

firms the generality of eq. fl36|) and provides strong support for the model 
developed in sect. [TTl Let us emphasize that in the present case the predic- 
tion for r2(a;) is obtained without any free parameter since (J2, ri, Rq, and 
e are known experimentally and is found by nonlinear rheology. 



C. Discussion and perspectives 

Our main result is that superposition rheology can be used to infer con- 
clusive information on the dynamics of wormlike micelles in the shear band- 
ing regime. In particular, superposition measurements lead to an estimate 
of the velocity cq which characterizes the dynamics of the interface between 
shear bands. The present uncertainty on the estimation of Cg through su- 
perposition rheology alone is of the order of ±40%. In our opinion this 
relatively large uncertainty is due to the use of a Mooney-Couette cell and 
to the subsequent approximations needed to process the superposition data 
in order to recover cq. Experiments in the cone-and-plate geometry under 
controlled shear rate should be simpler to process and should provide a bet- 
ter accuracy on cq. To minimize boundary effects in the concentric cylinder 
geometry, one could also avoid the use of a Mooney-Couette cell by trapping 
an air bubble below the inner cylinder. 
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FIG. 12: Amplitude of the interface displacement | r2 | versus lo. The solid 
line represents the prediction of eq. (j36p where all the parameters (T2 = 1 Pa, 
Uc = 68 Pa, ri = 0.5 mm, and cq = 0.28 mm.s~^ are known independently. The 
dotted lines were computed using cq = 0.25 (lower curve) and cq = 0.31 mm.s~^ 
(upper curve) in eq. (f36l) . They illustrate the sensitivity of the prediction to a 
10% variation in cq, which corresponds to the experimental uncertainty on the 
fit of fig. [TOl The fluid under study is a 6% wt. CPCl-NaSal solution. 

Let us now discuss the value of cq f ound from the superposition ex - 



periments reported above. According to [Radulescu and co-workers 2003 
Co is linked to the diffusion coefficient T) of the stress across the stream- 
lines, a central parameter in re cent theoretical approaches of shear banding 



Olmsted and co-workers 2000| . More precisely, one has 



da 
dc 




- = KGoyl^, (39) 
Co 



where c = eda/dt is the velocity of the interface, K is a. dimensionless pa- 
rameter that depends on the constitutive model. Go is the plateau modulus, 
and Ti the main relaxation time already introduced in sect. II B[ Following 
Radulescu and co-workers 2003j ]. we take KGo/riUO,'ji)'j]: = 0.3. From the 



nonlinear rheological measurements of sect. II Bl and using eq. (I7j), we find 
V\\{0,ii)ii ^ O.OQdc. With Co = 0.1 mm.s"^ and n = 0.87 s, eq. (ESI) yields 
V ^ 6.3 10~^2 ni2.s-i for the 8% wt. CPCl-NaSal solution. This corre- 
sponds to a stress correlation length ( = \J'Dt\ ~ 2.3 /xm, which is much 
larger than the mesh size ^ ~ {kT / GqY^'^ ~ 30 nm of our system. Let us 
emphasize the fact that ultrasonic velocimetry has provided the same order 
of magnitude for co at a slightly lower surfactant concentration but with a 
much better accuracy (co = 0.28±0.03 mm.s~^ for a 6% wt. CPCl-NaSal so- 
lution), which confirms that the stress diffusion coefficient estimated from 
eq. fl39l) should be in the range 10~^^-10~^^ m^.s~^ in the system under 
study. 
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Such a value of T) differs by two orders of magnitude from the stress 
diffusion coefficient inferred from transient rheo-optical measurements by 
Radulescu and co-workers 20031 ] in various wormhke micelle solutions with 



0.3 M CTAB {V ~ 1.2 - 7.2 10"^^ m^.s"^). Consequently our stress 
corr elation length is about twenty times larger than the estimate found 
by |Radulescu and co-workers 2003 |. C, — 100 nm, which was compara- 



ble to the mesh size of their micellar network ^ 26 nm). Since we 
used the CPCl-NaSal system rather than CTAB solutions such a difference 
may not be too unexpected. For instance the viscosity difference between 
the two coexisting phases, and hence the width of the stress plateau, is 
much smaller in our case (rf^/y — 7/) ~ 4.5) than in the experiments of 
Radulescu and co-workers 2003j | {t{'^n — ii) — 19-80). Moreover the value 



of K is not only rn odel-dependent but could also vary with the average shear 



rate [Dhont 1999[ |. Estimates of V inferred from eq. (l39i) should thus prob- 
ably be taken with care. In any case, superposition experiments in CTAB 
solutions, where wider stress plateaus and a better precision on their limits 
should make the determination of cq more accurate, would be very useful iri 
order to confirm the values of V found by [Radulescu and co-workers 20031 
in this system. 

Besides extracting a value for V, superposition rheology in the shear 
banding regime would be even more interesting if it could provide some 
information on the rheological behaviour of the nematic phase as suggested 
in sect. Ill A[ Indeed information about the structure and dynamics of the 
shear-induced, oriented phase is often tricky to derive from conventional 
measurements due to the slope in the flow curve that results from curva- 
ture and due to instabilities that occur on the high-shear branch. Here the 
dynamical behavior ?7j|^(u;) of the shear-induced phase may be recovered 
by considering the slopes of the linear fits of l/r^y in fig. [5] which are equal 
to i/vli ~ l/''7||Ar according to eq. (j33l) . The reconstructed ?7j|^(a;) data are 
presented in fig. Mjo), where they are compared to superposition data mea- 
sured at the beginning of the high-shear branch of the flow curve. The fact 
that the experimental data are systematically lower than the reconstructed 
data can be easily explained by the distance from the experimental shear 
rate (7 = 11.1 s^"*^) to the upper limit of the stress plateau {•jiy ~ 7.4 s"""^). 
Moreover, as shown in fig. Et^a), which compares experimental data recorded 
just below 7/ and the r]*^j{uj) data reconstructed using the ai = limit and 

Co = 0.1 mm.s~^ in eq. ( l33l) . superposition measurements at the onset of 
shear banding also yield a very good approximation of the complex viscosity 
?7yj(a;) close to the beginning of the plateau. This allowed us to check the 
consistency of our fitting procedure and to confirm that superposition rhe- 
ology provides useful quantitative information on the dynamical behaviours 
of both the entangled and the oriented states. A deeper analysis and mod- 
elling of such behaviours are left for future work. 

For the sake of completeness, fig. [13] shows superposition data obtained 
on the high-shear branch of the flow curve. Although a simple interpretation 
of fig. [T3r b) may not be possible due to the occurrence of flow instabilities 
for 7 ^ 15 s~^, these data clearly show that the dynamical behaviour of the 
shear-induced phase totally differs from the initial Oldroyd-B behaviour of 
the weakly oriented, entangled phase. Such information may turn out to be 
crucial for the modelling of shear banding since the exact behaviour of the 
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FIG. 13: Superposition rheology of an 8% wt. CPCl-NaSal solution in the high- 
shear regime: | r7jj(w,7i) | versus uj for (a) 71 = 11.1 and (b) 15.0 s~^. The sohd 
hnes correspond to an Oldroyd-B fluid (eq. (fTOl) ) with rjQ = 122 Pa, si = 0.59 s, 
S2 = 0.13 s, Ti = 0.87 s, and T2 = 0.60 ms. 

fluid at the limits of the stress plateau is usually unknown. 

Finally the influence of normal stresses or flow-concentration coupling 
in superposition experiments and the way to include them in a model also 
constitute directions for further research. 



Conclusion 

In this paper we have shown that superposition rheology constitutes a 
useful tool to access the dynamics of a shear-banded flow. A two-fluid semi- 
phenomenological model was proposed based on the simplest shear banding 
scenario. This model was shown to provide a good description of the oscilla- 
tions of the interface between shear bands in CPCl-NaSal wormlike micelle 
solutions sheared in the Mooney-Couette geometry. In particular an esti- 
mate for the stress diffusion coefficient T) was reported for the first time in 
the CPCl-NaSal micellar system, whose value was shown to be significantly 
larger than that reported for CTAB systems. Independent measurements 
of the interface dynamics through local velocimetry experiments nicely cor- 
roborated our model without any free parameter. We have shown, however, 
that a more accurate determination of the characteristic velocity requires 
a simpler and better controlled geometry. Further experiments, e.g. under 
controlled shear rate in the cone-and-plate geometry, should allow one to 
probe even more precisely the dynamics of the shear bands using only a 
standard rheometer and to infer important information on the dynamical 
behaviours of the two coexisting phases. The formalisms to use for these 
experiments are also supplied in the present work. 
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APPENDIX A: TWO-FLUID CALCULATIONS IN EXPERI- 
MENTAL GEOMETRIES 

In this appendix the detailed calculations for two-fluid superposition rhe- 
ology are presented in the standard geometries used in the experiments 
namely cone-and-plate, Couette, and Mooney-Couette geometries. 



1. Cone-and-plate geometry 

Let us first consider a cone-and-plate geometry of angle f3 <^ 1 and 
maximum radius Rq. In such a geometry and in a homogeneous fluid, the 
shear rate can be considered as constant throughout the sample. In the 
shear banding regime, 71 and ai are still linked by the lever rule (fT9l) so 
that eq. (120!) remains valid. Moreover, at a given distance r from the axis 
of the cone, the system is equivalent to an infinite parallel plate geometry 
of gap e = r tan /3 for which the shear stress is 

a(r)=a, + r/||(r)72e*-*, (Al) 

where ri*^{r) is computed from eq. (IMl) by setting e = rtan/5. One can then 
calculate the total stress exerted on the cone from 



Ro 

a{r)rdr. (A2) 



0^0 



In analogy with eqs. (125!) and (126!) let us define the two characteristic vis- 
cosities 

' '""^+^, (A3) 
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(A4) 



so that iri£) corresponds to the dynamical term 77/)^ of the complex viscosity 
(eq. with e = i?otan/3. With these notations, inserting eq. (lAip into 

eq. (1A2I1 leads to a = ac + ?7j| 72 e*"^* where 



Vw = Vl 



2^ + 24ln + ^ 
Vd Vd V Vl 



(A5) 
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The above expression for ?7y in the cone-and-plate geometry clearly differs 
from eq. (1241) obtained for infinite parallel plates. In particular a linear fit 
of l/r^y vs «! does not seem relevant. As discussed in sect. Ill A[ one could 

still use eq. f lASP to fit rj*^ with two free parameters 77^^(0;) and ijd^uj). Such 

a procedure would provide an estimate for rju^ui) and therefore Cq. Another 
way to proceed is to notice that in our experiments rjo > tjl, so that 

112 , 
-^- + — , (A6) 
^7^ Vl Vd 



which is equivalent to eq. ( 1241) with 2e = Rq tan l3. In this case the linear 
regression of 1/?7m may also lead to a good approximation of Cq. 



2. Couette geometry 

Let us now consider a concentric cylinder geometry (Couette geometry) 
where the inner cylinder of radius Rq is rotating while the outer cylinder 
of radius -Ri remains fixed. This choice is made to be consistent with the 
experimental section but our model can easily be adapted to any rotational 
configuration of the two cylinders. The gap between the rotor and the stator 
is e = i?i — -Ro- In the Couette geometry the shear stress is not homogeneous 
throughout the whole cell. Under the steady-state approximation already 
discussed in sect. Ill Al the shear stress depends on the distance r from the 
inner cylinder as 

^(^) = 7T ' 

1 + ^ 1 



Ro J 

so that cTi + (T2e"^* = cr(0) corresponds to the shear stress at the inner 
cylinder. Note that the rheometer may rather indicate "average" shear 
stresses for o"i and (T2 measured in superposition experiments. Since these 
stresses only differ from the values at the inner cylinder by a geometrical 
factor of order 1 and since this factor also depends on the way the average 
is defined, we shall leave out this complication and stick with ai and a2 as 
the values at the inner wall. 

In the simple shear banding scenario described in the introduction, the 
inhomogeneity of a ensures that only two shear bands separated by a single 
interface coexist in the gap. More precisely the shear-induced transition oc- 
curs when there exists < < e such that cr(rc) = (Tc- For r > rc, 0"(r) < cXc 
so the fluid remains entangled and in the high-viscosity state, while for 
r < rc, cr(r) > and the fluid is in the shear-induced low-viscosity state. 
Another c onsequence of the stress inhomogeneity is that the stress pla teau 
is not fiat [Radulescu and Olmsted 2000l . ISalmon and co-workers 2003 |. In- 



deed the shear-induced state first appears when ct(0) = ac and fills the whole 
cell when a{0) = adl+e/RoY- When e/Ro <^ 1, this leads to a linear a vs 7 
curve with slope da/d'y = 2e(Tc/-Ro(7/""7Af)- In the case of the experimental 
data shown in fig. [21 c/Rq ^ 0.04 and the shear stress is indeed seen to in- 
crease linearly in the shear banding regime. However, as already noted, the 
high-shear branch of the flow curve is hardly distinguishable from the "stress 
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plateau." Still we can take advantage of the existence of a tilted plateau to 
estimate 'Jn- Fitting the flow c urve at high shear rates by a B ingham fluid 
a = as + rjsj (as suggested by (Salmon and co-workers 2003| ) and looking 
for the shear rate corresponding to a = ac = 100 Pa yields 'Jn = 7.4±0.4 s~^ 
(see dashed line in fig. [2]). 

Thus, from eq. (1A7I) . it is required that Cc < cri ± 0-2 < 0-^(1 + e/ RqY 
for a superposition experiment to be performed in the shear banding 
regime at all times. Let us define ri such that = cri/(l + ti/RqY 
and Trlt) the position of the int erface at time t. The model proposed by 
Radulescu and co-workers 19991 ] implies that 



Co dt (Tc V '^1 / V -^0 + r. 



2 



Idr. a(r,)-a. ^ + ^ ) ( 1^ ) - 1 • (A8) 



In the linear response, eq. (lASP leads to rc(t) = ri + r2exp(za;t), where 



Ro+ri 

Since rdt) = a(t)e, one gets a(t) = ai + a2exp{iut) with 

= = (AlO) 

e e VV '^c / 

^2 ^0 Co 

0^2 = = -, TTT n . All 

e a, (/2o + ri)2 zcue + ^ ^ ^ 

This last equation is tested experimentally through velocity profile mea- 
surements in sect. IIII B[ 

Once the interface motion is known from eqs. (lAlOP and (lAllI) . one can 



go back to the apparent shear rate, i.e. the shear rate averaged over the 
whole sample 

71^) = 7i + 726 - ' 



\^^(^i('^)) V\\Ni^,(^i{r)) ^ 

+ r (A12) 

where o"i(r) = o"i/(l + r / RqY and cr2{r) = 0"2/(l + r / RqY . Since our super- 
position experiments are performed under controlled stress, we have noted 
the viscosities r//, 777V, 77^, and r/jj^ as functions of the local steady shear 

stress cri(r). In principle, knowing the different viscosities (from experimen- 
tal measurements or extrapolated data as mentioned in sect. Ill A|) . eqs. f lA9P 



and flA12p allow one to solve for 71 and 72 and thus to find ?7i| = (72/12- 



In order to get an explicit form for 77^ that we may compare to eq. 
we shall assume that the small-gap approximation e <^ Rq holds, which is 
almost always the case in standard experiments in the Couette geometry. 
In that case, expanding eq. ( 1A12I) to first-order in e/Ro and looking for the 
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constant terms leads to 



1 _ 
V 

^ e 

Rq 



7i 
0-1 



+ 



(A13) 



rii{(yc 



■ da '"^ 



l + af 



N I 



VN[(yc) 



da 



+ 1 



This yields the apparent viscosity t] indicated by the rheometer in the shear 
banding regime (up to some multiplicative factor of order 1 that depends on 
whether the rheometer actually indicates the shear stress at the inner wall 
or some average shear stress, as already mentioned above). Note the first 
order correction in e/Ro to the case of simple shear given by eq. (l20l) . By 
looking for the terms proportional to exp{iujt) in the first-order expansion 
of eq. (1A12I1 , one finds 



1 



72 

cr2 



1 1 1 
- + — + — , 

Vd Vd 



(A14) 
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where we have dropped the dependence on ac of the various viscosities for 
the sake of clarity. Equation flA14p generalizes eq. to the case of a 
small-gap Couette geometry and shows that ?7j| now involves three terms: 
the lever rule rji, the dynamical component rjD that arises from the motion 
of the interface, and rjQ a first-order correction to rji similar to that found in 
eq. (IA13P and linked to the stress inhomogeneity. Keeping in mind that the 
various viscosities in eq. ( 12^ are taken at ai = ac, the case of two infinite 
parallel plates is easily recovered from eqs. flA13p - flA17p when Rq oo. 

Therefore, in a small-gap Couette geometry, l/?7i| is a second-order poly- 
nomial in «! whose value for ai ^ is 



1 1 
lim — = h 

V*\ V\\i 



In - li 



Co 



luje 



2coe 



(A18) 



This is very similar to eq. (1271) so that the same data analysis should lead 
to the measurement of fjoi^^) and to an experimental determination of cq. 



3. Mooney-Couette geometry 

Experimentally, in order to minimize boundary effects due to the finite 
height of the cylinders, one often uses a composite geometry, called the 
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Mooney-Couette geometry, made of a Couette cell of gap e with a cone- 
shaped bottom such that e = i?otan/3 ^ RoP- In a Newtonian fluid and 
in the small-gap approximation, this geometry ensures that the shear rate 
remains constant over the whole sample. 

Using the results obtained in sect. lAll and IA2t one can easily construct 
a model for superposition experiments in the Mooney-Couette geometry of 
height h by considering the proportions eco = (l + -Ro/2/i)"^ and ecp = 1 — 
of the surface respectively covered by the Couette (co) and by the cone-and- 
plate [cp) geometries relative to the total surface. The total shear stress is 
then simply given by a = ecoO'co + ^cp o'cp, which yields 

V = ecor]co + ^cpVcp, (A19) 
V\\ = (^coV\\co + (^cpV\\cp, (A20) 

where ricp and 77*^^ are given by eqs. (120|) and (]A6p . and ?7co and r]*^^ by 

eqs. (1A13I) and ( 1A14I) . To close this set of equations, one has to specify 
the values of ai in the two parts of the geometry. Since the shear rate is 
perfectly homogeneous in the cone-and-plate, the steady component of the 
shear stress acting on the cone is aicp = <Jc so that the steady component of 
the shear stress acting on the inner cylinder is aico = {cri — ^cp o"c) / Cco- Thus 
the local proportions of shear-induced structure aicp and aico are given by 



«ico = — U/ ^ , (A21) 

e VV ecoO-c / 

aicp = . (A22) 

In - li 

In the limit c/Rq <^ 1 one can define an effective ai for the whole cell: 



Ro ( (^i-^cpCTc ,\ , 71-7/ /AoQ^ 

tti = ecoaico + ecp aicp = ^co — \ I ] + ecp — ■ [A26) 

e V V ecoac J In - li 



With eqs. flMOl) . fM . flMil) . flMT]) . and flX22l) . one can in principle 



determine the characteristic velocity cq and the dynamical behaviours of 
the two coexisting phases 77^^- and ?7y^ by fitting 77^ using eq. flA20p at a 
fixed oj. However, in practice, such a fit requires to know precisely cTc, 77, 
and 7Ar together with 71 and 77* for at least four different values of ui. As 
already pointed out, cxc, 7/, and 77V may be difficult to access and, in a curved 
geometry, are known to within 10% at best. Therefore the complexity of 
the fitting procedure along with the high number of unknowns prevent us 
to fit experimental data to the full model described above. Moreover the 
simple data analysis proposed in sect. Ill AI and based on an extrapolation 
to ai = (in order to remove the dependence on the unknown viscosity 
ryjl^) is no longer possible in the Mooney-Couette geometry since aicp and 
aico do not go to zero for the same 71 or ai. Nevertheless, in sect. Ill Bl it is 
shown that eq. (|24l) along with ai calculated from eq. (1A23I) may still allow 
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us to estimate cq. 
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